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Abstract. We prove that the Birkhoff sum S n (a) = — Ysk=i g(ka) of g = 
cot(7ra;) and with golden ratio a converges in the sense that the sequence of 
functions s n (x) = S[ X q 2n \/q2n with Fibonacci q n converges to a self similar 
limiting function s(x) on [0, 1] which can be computed analytically. While 
for any continuous function g, the Birkhoff limiting function is s(x) = Mx 
by Birkhoff 's ergodic theorem, we get so examples of random variables X n , 
where the limiting function of Sy xn -\/n — > s(x) exists along subsequences for 
one initial point and is nontrivial. Hardy and Littlewood have studied the 
Birkhoff sum for g'(x) = —ttcsc 2 (itx) and shown that -\ Y^k"i 9'(^ a ) s t a y s 
bounded. Sinai and Ulcigrai have found a limiting distribution for g(x) if both 
the rotation number a and the initial point 9 are integrated over. We fix the 
golden ratio a, start with fixed 8 = and show that the rescaled random walk 
converges along subsequences. 

1. Introduction 

We study the sum S n — J2k=i ff(^ a ) with unbounded g(x) and irrational a. If 
the function g is not in L^T 1 ), there is no ergodic theorem to compute the limit 
of the Birkhoff averages S n (a)/n. We will see however that there is a convergence 
in a rather particular way if a is the golden mean and g has a single pole like for 
g(x) = cot(irx). We have studied this particular Birkhoff sum for the antiderivative 
G(x) = log \l-e 27Tlx \± of 5 in [TJ with the help of Birkhoff functions s m (x) = %^ 
on [0,1]. The function s m (x) encodes the random walk {SkYk=i- For identically 
distributed L 1 random variables, Birkhoff 's ergodic theorem Xk(6) = g(9 + ka) 
assures that s m (x) — > Mx, pointwise for almost all 9, where M = E[X] = L g(9) d9 
is the expectation. We call the limit s(x) = Mx the Birkhoff limiting function 
of g- 

Our main result is to that that for the random variables Xk(9) = cot(7r(6' + ka)) 
and particular starting point 9 = 0, the functions s rn (x) converge pointwise to 
a selfsimilar Birkhoff limiting function s(x) along subsequences q n , where p n /q n 
denote partial fractions of a so that q n = p n +i are the Fibonacci numbers: 

Theorem 1.1. The Birkhoff limiting function 



s(x) = lim 
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for g(x) = cot{~Kx) exists pointwise along odd and even subsequences. The graph of 
s(x) = lim n _>oo S2n(x) satisfies s(ax) = —as(x) and is continuous from the right. 

We call the graph of s the golden graph. It can be seen in Figured! We will de- 
scribe the function s(x) as an explicit series which can be evaluated in 0(log(l/e)) 
steps if we want to know it up to accuracy e. An analytic function a will be key to 
determine s. 

Why is this interesting? Some motivation and references were already given in 
PHI [11], where we studied the Birkhoff sum for the function G satisfying G' = g. 

The Sinai-Ulcigrai model 

g(x) = 2(1 - exp(^))- 1 = 1 + i Sm ^ x ) =l + i cot(^) 

1 — cos(x) 2 

is equivalent. In |13) these authors have shown the existence of a limiting distribu- 
tion of S n /n if one averages both over the initial point 9 and the rotation number 
a. In other words, they have shown that the Birkhoff sum for cot and the map 
T[0, a) = (9 + a, a) on T 2 has a limiting distribution in the sense that the random 
variables S n /n converge in law to a nontrivial distribution. Where Sinai-Ulcigrai 
average over all initial conditions and rotation numbers, we look at one specific 
orbit and take one specific rotation number a, the golden ratio. We prove that in 
this particular case, the rescaled random walk S q2rl /q2 n converges. For the Sinai- 
Ulcigrai distribution, it is not required to go to a subsequence, but one has to 
integrate over 9 and a. 

The related aspect of the Sinai-Ulcigrai limiting theorem will be touched upon 
at the end of the paper, when we look at the antiderivative G of g, the Hilbert 
transform of the piecewise linear Hecke function H studied by number theorists. 
Since Denjoy-Koksma theory 015] implies that the log Birkhoff sum of S„ /log(n) is 
bounded for almost all 9 and a of constant type, also the Birkhoff sum S„ / log(n) for 
the Hilbert transform G has this property and there is a subsequence S^ k / log(rife) 
with a limiting distribution. We have explored this limiting distribution for a sin- 
gle orbit in in [11] . The result proven in the present paper implies that for the 
golden ratio a, the Birkhoff sum S% n /q k has a limiting distribution for derivatives 
g( k \0 < k, again for a single orbit. 

In probability theory, the Birkhoff sum under consideration deals with examples 
of ergodic aperiodic sequences of identically distributed random variables X^ with 
finite Cauchy principal value expectation. We can so give explicit expressions for 
the sum of Cauchy distributed random variables 

S n {6) = X 1 {0) + ---+X n {6) 

for one point 9 — 0. This is related to Sinai-Ulcigrai, who take as the prob- 
ability space the two torus and generate the random variables with the shear 
T{9, a) = {9+a mod 1, a) and find a limiting distribution when averaging over pairs 
(9, a) of initial conditions and rotation numbers. Because every natural probability 
space is equivalent to the unit interval 57 = [0, 1] with Lebesgue measure P = dx, 
every sequence of identically distributed Cauchy random variables can be obtained 
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Figure 1 . The golden graph is the graph of the Birkhoff lim- 
iting function 

s(x) = lim S [q2nX] /q 2n = - lim S [q x] /q 2n +i . 

n— too 1 1 n— >oo L J 

This limit is more interesting than for continuous functions g, 
where by Birkhoff's ergodic theorem, the Birkhoff limiting func- 
tion is a linear function s(x) — Mx. 



as Xk = cot(7rT fe (a;)) with some measure preserving transformation T : Q — > f2. 
Taking an irrational rotation gives stochastic processes with strong correlations be- 
tween individual random variables. The "risk" of this process is large because the 
orbit gets close to the pole, where the "loss" or "gain" can become arbitrary large. 
It complements a trivial "Cauchy central limit theorem" which tells that for any 
process X n = Y n + Z n with identically distributed but not necessarily uncorrelated 
L 1 random variables and IID distributed Cauchy random variables Zk, the sum 
S n /n converges in law to a Cauchy distributed random variable. This is 'triv- 
ial' because (1/n) Y^k=i ^ k conver g es to a constant M and because (1/n) Y^k=i 
all have the same Cauchy distribution. The Cauchy distribution shares with the 
Gaussian distribution the stability of adding independent random variables. Since 
the standard deviation is infinite, the Cauchy distribution can model high risk. In 
this paper, we study a limit theorem in the case when the Cauchy process has corre- 
lations given by a Diophantine rotation. Since the result will be pretty independent 
of the function g as long as it has a single pole and the initial point is that pole 
it can - like the Sinai-Ulcigrai result - be seen as a limit result, where a chaotic 
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Bernoulli process is replaced by a Diophantine irrational rotation. 

Theorem (jl.ll) allows to estimate S n /n as (1/n) 53/=i s ( x i)i where the Zeck- 
endorff representation n = Y2k=i a k1k of the integer n is linked to the (3- 
expansion xi = X)fe=i a k<^ L k - Similarly than in the bounded case, we could 
give an estimate of S^q 100 f° r example, even so there are too many terms to add 
this up directly. The reason is that we can compute S n for Fibonacci lengths 
as S n ~ ns(n/q), where q is the smallest Fibonacci number larger or equal to n. 
Given a number n like n = 23'866, we can represent it as a Zeckendorf sum 23866 — 
17711 + 4181 + 1597 + 377 of of Fibonacci numbers (722,919,917,914- The next larger 
Fibonacci number is 46368 = 924. The real number x = 23866/46368 = 0.514708 
has the f3 expansion x = a 2 + a 5 + a 7 + a 10 , where a is the golden ratio. We 
will see that S n ~ X)fc°=o a k(~ a ) k<J {yk) with explicitly known y^. In our example, 
V\ = 0,2/2 = -0.0344396,1/3 = -0.026314; j/ 4 = 0.00118918. While the actual sum 
is —0.187542, the analytically computed sum obtained by adding up 4 summands 
is —0.18763. We have made use of a function a which n independent and analytic 
and which can be computed at first. While S'xgioo/lO 100 can not be computed by 
summation, we could determine the sum analytically. The accuracy depends on the 
structure of n, on how many Fibonacci numbers are needed to Zeckendorff represent 
the number. As the function a is analytic, the Birkhoff sum can so be estimated 
so fast that it could be called "integrable" and be diametrally opposite to the IID 
case, where the value of S n /n can only be computed in average for numbers like 
10 100 . 

For number theory, the Birkhoff sums under consideration are related to the 
values of a modular form at a specific point. The sums S n always converge without 
normalization if a is in the upper half plane {a G C | im(a) > }. The classical 9 
function 

00 

o(a) = e "™ 2Q = w " 2 = 1 + 2 J2 w " 2 

neZ n£Z n=l 

with w = e™ a is related to the modular form 

00 

R{a) = 9 2 (a) = l + 2^c7(fca)+3(/ca+7r/2),= l+2S 00 {co\, i 7r/4)-2S' 00 (cot, — tt/4) 
fc=i 

for im(a) > 0. In other words, the modular form R is up to a constant the sum of 
two Birkhoff sums described in this article which converge for a in the upper half 
plane. 

The sums also relate to the theory of partitions in additive number theory 
because g(x) = G'(x) with G = log(l — z) and z — e 27rlx . The Birkhoff sum for G 
is the logarithm of a product 

fe=l 

where w — e 2ma . The function — G is related to the generating functions for the 
partition function p(n) because P^ 1 = J2k=iP(^) wk + 0(w n+1 ). 
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To illustrate the connection between number theory and dynamical systems, lets 
prove a well known statement in the theory of integer partitions p(n) (which tells 
in how many ways we can write n as a sum of positive integers) with the help of 
Birkhoff sums in dynamical systems theory: 

The partition function p(n) satisfies limsup„ ^(w)! 1 ^™ = 1. 

Proof. Euler's pentagonal number theorem (i.e. |4j) shows that the function 
Q( z ) — IlfcLi(l — zk ) = Sn a ™ z " which is related to the partition numbers by 
P(z) = J2n=iP( n ) zn = ft^LA 1 ~ = 1 /Q( Z ) has Hadamard gaps. By 

Hadamard's lacunary series theorem, the Taylor series has a natural bound- 
ary as a circle. This implies that establishing boundedness on one point \z\ = r 
implies that the radius of convergence is larger than r. Define g(a) — (1 — re 27 ™") -1 
so that /(a) = log(P(z)) = — J2T=i 9^°) For r = \z\ < 1, and z — exp(27ria) with 
Diophantine a, the theorem of Gottschalk-Hedlund [9 shows that the function / is 
a coboundary f(x) = F(x + a) — F(x) implying that the sum is bounded. Having 
so proven that the Taylor series has no singularity on \z\ < l,the partition function 
p(n) satisfies limsup n |p(ti) | x / rl > 1. The other inequality limsup^ ^(n)! 1 /™ < 1 is 
trivial because p(n) > 1 implies Ysn=i p{ n )z n > Yln=i z ™- ^ 

While this growth-statement for the partition numbers p(n) is not difficult to 
prove directly |12j . the just given new argument illustrates how complex analysis, 
dynamical systems theory and Diophantine notions can play together. 

The choice of the core function g (x) = cot(7ra;) for the Birkhoff sum is distin- 
guished because the periodic function g has as a distribution the constant Fourier 
expansion 

oo 

(1) g(x) = cot(7rir) = 2 sin(27rfcx) 

fc=i 

as can be seen by taking the real part of 

ix , -ix i I 2ix °° 

cot(x) = i . + . - -i^V = + 2 T e 2lkx ) . 

k=l 

How special the cot function is can also be seen with the Euler's cot formula 

(2) sw = -i;^t. 

7T * — ' x — k 

Appearing also as the Hilbert kernel it is the most natural periodic function with 
a single pole. 

In [TDJ[TT], we studied the Birkhoff sum of the antiderivative of g(x) = cot(7rx), 
which is 

G{x) = log(2-2cos(27ra))-!- 

2ir 

= log (2 sin(7rx)) — 

7T 

= log|l-e 27rix |- . 
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Figure 2. The rational golden graph is the Birkhoff limiting 
function 

*(*)= lim ^I(T n ) 

in which the translation T n {x) = x + P2n/q2n is periodic. It is 
different from the Birkhoff limiting function s in which T(x) = 
x + a is n independent. The function t is not afhne self similar 
(a + b(—t(ax) — t(a(l — x)) only comes close). But the function t 
can be useful to answer questions posed in [TT] . 



Euler's formula for g(x) = G'(x) can also be deduced by logarithmic differentiation 
of the Euler's product formula for the sinc-function 

, . sin(7rx) t-t x — k 
smc(7ra;) = = — - — 

X7T k 

fc#o 

because 2sin(7rx) = exp(7rG(a;)). 

The current paper will lead to explanations to some of the statements found in 
the experiments done in [11]. Wee do not give much details here yet but say that 
one connection is the fact that the difference of the Birkhoff sums of the antideriva- 
tive G(x) of a; for successive continued fractions p n /qn,Pn+i/qn+i can be expressed 
by Rolle's theorem as a Birkhoff sum of g' but with a different initial condition. 
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While the golden graph of g is selfsimilar, the limiting distributions seen in [II] are 
only almost self similar but the difference can be expressed by a Birkhoff sum for 
g' which is monotone. 

While the rotation number a is special because it has a constant continued frac- 
tion expansion, the story could be generalized a bit. For quadratic irrationals, there 
is a periodic attractor for the limit (see pTj). 

For L 1 functions /, by the Birkhoff 's ergodic theorem, the limiting function 
is for almost all 9 given by 

[xn] m 

s(x)= lim - ^Tf(T k 8)= lim — V f(T k 9) = Mx , 

fe=l fe=l 

where M = j f(9) d9. We are not aware of any other case of a stochastic process, 
where the Birkhoff average s n (x) = S\ nx ]/n converges to a nonlinear function on 
[0, 1] along an orbit. The cot function can be modified by adding any continuous 
periodic function F: 

Corollary 1.2 (A single orbit ergodic theorem for the golden rotation). For any 
function g on T 1 = [0, 1) = R/Z which is continuous away from 9 = and has a 
single pole at 9 = 0, the Birkhoff limiting function at 9 = 

8 ( X ) = lim 

n->oo q 2n 

is the sum of a multiple of the golden graph c ■ s(x) and a linear function Mx. 

Proof. Because / has a single pole at 0, we can write f(x) = ccot(7ra;) + F(x), 
where F(x) is continuous and especially in L . The Birkhoff limiting function for / 
is the sum of the golden graph cs(x) and the Birkhoff limiting function of F, which 
is Mx with M = F(x) dx. □ 

We will be able to compute the function six) explicitly from a function 




which we will show to be analytic in y. Its Taylor expansion starts with 
o-O) = 0.258 + -1.2472% + 0.6736?/ 2 - 1.38082y 3 + ... . 
First, we look at the analogue function for rational a, where the function is 




The function r is odd and will be given explicitcly: 




From a we can then get s. 

An indication why things work nicely is that the cot function and their deriva- 
tives have structure. Hardy and Littlewood have made use of this early on. The 
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Figure 3. A caricature for the golden graph is the function 

oo 



fe=i 



where a = 1) /2 and a k is the /3-expansion of x. The function 

o has the same symmetry o(a n x) — (—a) n o(x) as the golden graph. 
The golden graph is a modification of this graph since we will show 

oo 

s i x ) = ^a k {-a) k a{y k ) , 

k=l 

where y = 0, y fc = lim n ^. 0O [|a;fcg T ,|a]gr n with x fc = J2i=i a i al - 



selfsimilar nature of the Birkhoff sum appears in the following statement which we 
have used already in integrated form in [10] and which will be used in an equivalent 
way below in Lemma [ 



Proposition 1.3 (Birkhoff renormalization) . For any periodic translation T p / q (jj) = 
y + p/q on the circle R/Z with gcd(p, q) = 1, the function g(y) = cot(7ry) is a fixed 
point B p / q (g) — g of the linear Birkhoff renormalization operator 



Remark. 

We believe that any periodic function which is a fixed point of all the Markov 
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Figure 4. The Markov operator B 1 / z is illustrated in the hgure 
as it is applied to a general function, where it produces a new 
different function. The two figures below show the operator applied 
to the cot(7ra;) function which is a fixed point: 

1 y y 1 y 2 

-[C0t(7r-) + COt(7T- + 7T-) + COt(7T- + 7T-)] = COt(7TJ/) . 

We see the three graphs, obtained by taking stretched pieces of 
the cot graph and then averaging them to get back the function. 
Each operator B p / q is a Perron-Frobenius operator for an interval 
exchange transformation. 



operators B p / q must be of the form f(y) = a + bcot^iry) with some constants a, b. 
But we have not yet been able to prove such a fixed point result. 

Proof. The operators B p / q are Markov operators because they satisfy Bl = 1 and 
11-^/1 1 1 — I l/l 1 1 f° r positive L 1 functions. Use the cyclotomic formula 



-e 2mk / n z) = 




k=l 
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and plugging in z — exp(27riy), then differentiate with respect to y to get 

5-1 

Q 



1 \— * , (k + y) 

- > cot(7T ) = cot(iry) 



fe=0 



This is equivalent because {pk mod q } = {0, . . . , q — 1 } follows from gcd(p, q) = 1. 
We have so shown that g = a + 6cot(7ra;) is a fixed point of B p / q for all p/g. □ 



Remarks. 

1) Unlike for other parts of this article which often were found experimentally, 
the fixed point property was obtained first in a deductive way from the cyclotomic 
identity. 

2) Because related identities appear in formula collections, we must assume Propo- 
sition QT2] to be "known". Still, it is a remarkable gem and probably not "well 
known". Here are some related identities (1/g 2 ) X)fe=o csc 2 (7r(/c + y)/q) = csc 2 (y) 
or (1/g 2 ) X)fe=o c °t 2 ( 7r (fc + J/)/ ( 7) = cot 2 (?/) + l — l/q for cot 2 which are equivalent by 
differentiation and using the trig identity cot 2 (a;) + l = esc 2 (a;). We get so functions 
which are fixed points of the higher order Birkhoff operators 



1 fe=o q 



3) In the simplest case, the renormalization picture provides for p/q = 1/2 a double 
angle formula for the cot function: for x — ixyjl it is 

7T 

cot(a;) + cot(a; + — ) = 2cot(2a;) . 

For other p/q, we get "de Moivre type" trigonometric identities for the cot 

or tan functions analogues to the magic which de Moivre has provided for cos and 
sin. For example 

, , , 7T 2-7T. , 37T. . 47T 

cot(a;) + cot(a; H — J + cot(a; H ) + cot(.x H ) + cot(x H J = 5 cot(5a;) . 

5 5 5 5 

Before we launch into the details of the Birkhoff sum puzzle, here is a schematic 
overview over the proof. We start with an explicit Birkhoff sum in the rational 
case by using an elementary cyclotomic formula leading to the function r(y). In- 
terpolating this with Fibonacci steps allows us to get a related analytic function 
<y{y) for which we can compute the Taylor coefficients with arbitrary precision. 
This in turn provides us with a function s(x,y) via a j3 expansion of x. The slice 
s(x) — s(x, 0) is then the golden graph. This golden graph in turn allows us via 
Taylor expansion and Abel summation to compute analytic " golden graphs" in the 
rational case t(x,y). This function t is not selfsimilar but it will link us to the 
original observations done in [11]. It is interesting to note that the direct way in 
the rational case from r to t(x, y) is blocked because in the rational case, there is 
no self similarity. We need to make a detour over the golden ratio and the golden 
graph to understand the rational case! 
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Figure 5. The Zeckendorff representation of a number is the 
discrete analogue or dual to the /3 expansion of a real number. The 
picture illustrates the Zeckendorff representation ofl6 = 3+5+8 = 
Q2+Q3+Q4 as a sum of Fibonacci numbers as well as the j3 expansion 



of 0.763932... = a 2 
ratio. 



a + a as a sum of powers of the golden 



T(y)=t(l-,y) 



Steps 



t(x,y) 



Taylor 



t(x) = t(x,0) 



(Fig 2) 



c(y) = s(l ,y) 





r \ 

s(x,y) 


r 


) 


v. 



(Fig 7) 



/3-exp. (Fig 6) 



s(a;) = s(x, 0) 



Golden Graph (Fig 1) 



Acknowledgements: Thanks to John Lesieutre and Folkert Tangerman for 
previous joint work [KL] and [KTp] without which this paper would not have been 
possible. 



2. Constructing the function s 

Every real number x £ [0, 1] has a /3-expansion X^^Li a nOt n with ai £ {0, 1}. 
If a = (V5 — l)/2, then the possible sequences {a n } n are the ones for which 
no two consecutive are equal to 1 (see for example [3]). This is related to the 
Zeckendorf representation n = Yl^Lx a n1n of integers as a sum of Fibonacci 
numbers, where only finitely many a n are nonzero. Both in the /3 expansion of a 
real number in [0, 1] as well as the Zeckendorf representation of a positive integer, 
two successive l's are forbidden in the subshift of finite type of all {0, 1 } sequences 
because a k + a k+1 = a k ~ x and so q n -i + q n = Qn+i- Every point x e [0, 1) has 
a /3 expansion. Examples of some /3-rational points are x = a = 0.1, a 2 = 0.01, 



a + a fi — 0.101 Q or a 



a 



a 



0.10000101 Q . 



The semi- continuity of s is reflected by the fact that for every k > 0, the sum 
Sk+q n /q n has the same limit. For k < 0, the limit is different. To understand 
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Figure 6. The generalized golden graphs x — > s(x,y) are seen 
here for different initial points y. The figure shows the graphs for 
y € [—0.15,0.15]. We prove in this paper that the graphs deform 
continuously in y for small \y\ (one could actually verify it for 
\y\ < l/v5). With the initial conditions y n = y/q n with small \y\, 
the limit SLgj (y n ) — > s{x,y) still exists along even and odd sub- 
sequences. 



the limit, it is useful to start at different points and to look more generally at 
the Birkhoff sum S n (6) — Ylk=i d(@ + ka). We choose = yjqin and define the 
s-function 

s(x,y) = lim S [q2riX ](y / q 2n ) / q2n = - Urn S [q x] (y / q 2n +i) / q 2n +i ■ 

n— too n— too 

We show now that the function 

a(y) = s(l-,y) = Jim 5 [g2n+l] {^JL-J / q2n+l 

determines s(x, y) and that the graph of s is selfsimilar. We get 

Proposition 2.1 (The golden graph in terms of a). Assume we know that a(y) = 
s(l,y) exists. If x — Y^JLi ajOt? , then 



s(x,y) = ^2a k (-a) k a(y k ) , 



k=0 



where y = y, y k — y + lim n ^oo[\x k q n \a]q depend on x via x k = J2j=i a 3 a J . We 
have for every n the symmetry s(a n x,y) — (— l) n a n s(l, y). 



SELFSIMILARITY IN THE BIRKHOFF SUM OF THE COTANGENT FUNCTION 13 









0.4 


- 














0.3 
















0.2 
















0.1 










-0.15 


-0.10 


-0.05 






0.05 


0.10 


0.15 



Figure 7. The analytic function a(y) determines the function 
s(x,y) by a /3 expansion as shown in Proposition 12. 11 

Proof. Given x = cLiOt 1 define the real numbers Xf. = 2<=i a i& 1 ■ The selfsim- 

ilarity claim follows from s(ax,y) = —as(l,y) which is a consequence of the fact 
that for the golden ratio, 

±1 

aq n - q n +i = • 

q n 

Use that for x = a k , the selfsimilarlity implies that s(a k , y) = (—a>) k a(y(—a) k ). 
Use the notation [x] = x mod 1. Define yo — V and y^ = y + linin^oo |a;fc<? n | • a ■ q n . 
Now, 

oo oo 
fe=l k=l 

□ 

The main task is therefore to construct the function a. To do so, we have first to 
explore the rational case, where we sum up q n — 1 terms with the rotation number 
Pn/qn instead of with a. The analogue to the a function is called r. It can be 
computed explicitly for g(x) = cot(7rx) and defined as 

1 9-1 k 

This expression can be simplified and also explicitly be given in the limit q — > oo. 
Lemma 2.2 (The function r). If {kp }?Zi = {1, . . . , g — 1 } then 



T P /q(y) = g(y) - 
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In the limit q — > oo, 

r{y) = g(y) - — • 

ny 

Proof. The cyclotomic formula 

JJ(l_ e 2***/»3) = £_£ 

fc=l z 

holds for all complex numbers z =/= 1. Taking logarithms and looking at the real 
part, we obtain for 2 = e 2my and real nonzero y the formula 

£log|2 - 2cos(27r- + 2tt^)1 =21og| ? ~ | = G(y) 21og|e 2 ^ - 1| 

k—1 

so that 

n—l j 

£G(2 + £) = G(y)_G(J). 
' — ' n n n 

fe=l 

Now differentiate with respect to y to get 

n * — ' n n n 

fc=i 

Now use g{y/q) ~ q/(yn) to get the statement in the lemma. □ 
Remark. Repeated differentiation gives more identities 

n-1 



*T 4tE5 (0 (- + -)^/ ) ( 2 /)-/ ) (-)A 

n-s-oo rt l+1 ' — ' n n n 



in 1 

~n" 

fe=i 



For later, we also will need the rational analogue of the s function which we call 
the t-function 

[<k] 



t( x >y)= lim ~ y^g(~ + fc ~) 



We see its graph in Figure (0) for y = 0. We know that t(l,y) = r(y) is analytic, 
as we have given it explicitly. The function t is only close to self similar but does 
not quite fit an affine similarity t(ax, ry = t{x, y)a + b. Therefore, we can not get 
t from t as we did get s from a. In the rational case, we have an explicit limiting 
function for r x — l - but no self similarity. In the irrational case, we have a self 
similar s function but no explicit expressions for a. 

Corollary 2.3 (The i-function). If s(x,y) is continuous to the right and analytic 
in y then the function t(x,y) is continuous to the right in x and analytic in y. 

Proof. We can get t{x, y) from the function s(x, y) by making a Taylor expansion in 
a centered at the golden ratio and evaluate it at a continued fraction approximations 
p/q: 

t p / q (x, y) = s(x, y) + s'{x, y)£ - a)/l! + s"{x, y)(^ - a) 2 /2\ + ■■■ . 

q q 

and then take the limit p/q — > a. 
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Wehavc|H- a |<^and 

rl m 11 h m , , 11 

A m = ^g{ V - + ka)\ = — g ( m \V + ka) 
dy m \ n q m * \ ' 

We need to show that the Taylor series for t p / q (x, y) converges. For fixed q, it can be 
written as a sum J2 m A m /m\ which has terms which satisfy for large q (neglecting 
l/\/5<l) 

1 lqx] h m a 1 lqx] 

y fe=l y y y k=l 

With G; = X^=i 9 3 we can write the sum using Abel summation (discrete integra- 
tion by parts) as 

[qx] [qx] 

E ^ fe5fe = G[qx]h[qx\ ~ E ^feC 1 * - frfc-i) • 
fe=l fc=l 

Now 

1 ^ 1 I« x l frm l.m-1 1 

- E G fc (** - /*_0 = - E G fc (— - — ) = O(-) 

so that the second part after Abel summation goes to zero. With (1/ q m )G\ qx ] — > 
s( m )(x,y) and ^1 = (qx) m /q m = x m wc finally sec A m - s< m )(x,y)x m for g -> 
oo. " □ 

The following proposition gives a from the analytic function r: 

Proposition 2.4 (From r to <r). T/ie junction a(y) is analytic in y for small \y\ 
with Taylor expansion a(y) = Y^ilo a iV l / ^ > where 

1 9-1 

ai= lim -— Vfc 5 ( ' ) (- + M 

7l->oo g'+J- z — ' q 
k—1 

and q = q n is the n 'th Fibonacci number. 

Proof. We will compare the functions a to t which are given as the limits 

9-1 



1 

q k=i q 



and show that the difference is analytic. To do so, we write the difference as a sum 
of differences of higher order rational tau functions, then linearize 

o- n (y) - r n {y) = E T «+2'+2(y) - T n+2i(y) = E t' 2l+1 (y) ( 77^ - pp) + R n , 
7^o i=Z Ql + 2 Ql 

with 

> f \ 1 V" k 'tV , u Pl \ 

rnAy) = -T [ Q i 9(- + k- ) , 
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where q — q n ,P = Pn,Qi = qn+i,Pi = Pn+i for the sake of notational simplicity. 
The Taylor rest term 

R -VV 2) (v)(^- - —) 2 /2 

disappears in the limit n — > oo because ( J +2 —■§-) = (ttt=I — ) 2 ~ a 2 -^ ~ a 2+il -4- 
^ y Qi+2 Qi J y QiQi+2 J Qf 

so that remark after Lemma 12721 shows R n — > 0. Because T' n is a sum of negative 
terms, the boundedness of r' n ; (y) follows from the Cauchy-Schwarz inequality. We 
know that the limit n —¥ 00 of each r' n , {y) exists as an analytic function and that 
it is bounded above by Ca l . Therefore also the sum converges and is an analytic 
function. It is the difference between a and r we were looking for. □ 

Remarks. 

1) It is in this proof that we see the reason why we have to go along the subsequence 
<?2n to get the golden graph and not q n . The reason for the sign flip is: 

Pn+2 Pn +1 

= , n odd 

9n+2 qn 9n<7n+2 
Pn+2 Pn -1 

= , n even 

9n+2 9n 9n9n+2 

The t function has an analytic formula for any k, we do not even need a Fibonacci 
number. It is the difference between r and a which flips the sign. 

2) A Taylor expansion in a does not work when going from t to a. (We only make 
Taylor expansions with respect to y.) Abel summation indicates that this would 
require us to know Cesaro sums of the function t(x,y). The Taylor expansion will 
work however in order to get t from s. We have written therefore a — r as a sum 

f 1,1+2 °f smaller differences and estimated this as a Birkhoff sum with rotation 
number p n+ \/q n+ \. We can estimate each entry to be smaller than 1/qi and so 
get a Cauchy sequence. There is no need to use Cesaro convergence because the 
Birkhoff sums for g' contain only negative entries and establishing boundedness is 
enough. 

While we only need the first approximation, lets state a lemma generalizes 
approximations like \/5(a — p/q) = {— l) n+1 /q 2 + 0(q~ 4 ) and y/E(a — p/q) — 
(— l) n+1 /q 2 = l/(5g 4 ) + 0(q~ 6 ) etc. It is an exact formula estimating the error 
from the golden ratio a to any of its continued fraction approximations p n /q n '- 

Lemma 2.5 (Catalan gold). For any partial fraction p n /q n of a = (v5 — l)/2 we 
have 

r- V ^ ('-l % )("+ 1 )( fc + 1 )c l 1 

q n ^ 5 A qi k+l 

where Ck, — (2k)\ / (k\(k + 1)!) are the Catalan numbers. 

Remark. This lemma was found empirically by determining the coefficients 
Ck using high precision numerical methods using 1000 digit approximations of the 
algebraic numbers. Looking up the resulting integer sequence revealed the Catalan 
connection. Even so I know now an elegant proof, it might be well known, but no 
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0.8 - 



0.6 - 




Figure 8. The limiting function S[ xqn ]/q^ of the esc 2 Birkhoff 
sum is monotone because the summands are positive. A modifica- 
tion of this graph has appeared in as a discrepancy from the 
selfsimilarity of the difference graph. 



reference so far has emerged. For n = 1 for example, where q = l,p/q = 0, the 
Catalan gold lemma gives the formula 

The right hand side is c(l/5) = 2/(1 + - 4/5) = V5(\/5 - l)/2. The proof of 
the lemma has connections with Hilbert's 10th problem. The digging of Catalan 
Gold is only needed to find accurate bounds for the interval in which cr(y) converges 
and is left to the reader. 




3. Remarks and outlook 

We briefly look at the Birkhoff sum for the antiderivative G of g which was 
studied in [TT]. The function has relations with Hecke's investigations [7] because: 

Proposition 3.1. G(x) = log(2 — 2 cos(27ra))/2 is the Hilbert transform of H{x) — 
x- [x] - 1/2. 
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Proof. Writing down some identities makes this clear: 



. , , . . ^ sin(27rnx) 
H(x) = tt(x - [x] - 1/2) = - 22 — " " 



arg(l - e 2 ™) , 



oo 



cos(27rnx) 



G(x) = log(2- 2cos(27ra;))/2 = log |2 sin(7ra;)| = - 



n 



log|l-e 2,reE | . 



Note that while G'(x) — 7rcot(7ra;) = irg(x), the derivative of H is only defined as 
a distribution. There is no corresponding Hilbert dual result therefore for the cot 



Remark. The first identity for H{x) is geometrically amusing: if a point r(t) 
moves on a circle with uniform speed, then the direction vector from r(t) to r(0) 
changes with uniform speed. Since this can be done from a second point also, it 
implies Thales result about angles in a triangle where the base segment is fixed. 

As a direct consequence of functional analytic properties of the Hilbert transform, 
we get a rather general Denjoy-Koksma type stability result. It confirms some 
numerical experiments in [101 111] : 

Corollary 3.2 (Log boundedness of the Birkhoff sum). If a is Diophantine of 
constant type, then there exists a constant C such that for almost all x, we have 



Proof. The Hilbert transform is translational invariant. Therefore, the Birkhoff sum 
of the conjugate functions remains the Hilbert transform of the original Birkhoff 
sum. Since the Hilbert transform is a bounded operator on I? , we have | !$«(<?) ct)\ 1 2 < 
C\\S n (h, ck) 1 12 - By Denjoy-Koksma applied to to function h which has bounded vari- 
ation, the later grows only logarithmic. The former grows logarithmically too. □ 

We still have to prove all the statements in [11] but only outline here. Lets denote 
by the Birkhoff sum for G, where the summation goes from 1 to q n — 1. We focus 
again on the situation, where a is the golden mean. The methods developed here for 
cot will allow to show that the limit S„ exists. The difference limit function is the 
sum of a self-similar golden graph with different initial condition and a correction 
which comes from Abel summation and which is close to the graph of the esc 
Birkhoff limit. 

Theorem 3.3. The Birkhoff difference functions f{x) = lim n _ J . O0 S[ xqn ](p n +i/ 'q n +i)— 
Slxq^iPn/qn) of g{x) = log(2 — 2 cos(27t:e)) converge pointwise. For the Riemann 
sum of g{x) — log(2 — 2cos(ir)) 7 the limit S Qn converges and \og(S qn )/ log(q„) 
stays bounded. For the imaginary case, f2n{x) —> cx 2 ,fi n -\(x) — > —ex 1 with 
c = (1 + a 2 )7r/v5- Consequently, for the combination of real and imaginary part 
g(x) — 21og(l — e lx ) the difference limit exists. The limiting function f has the 
property that f{a) + a 2 f is monotone. 

Proof. The difference by a Rolle estimation, which is a Riemann sum for the de- 
rivative cot, but with a different starting point and a correction. To estimate the 
difference function, we estimate the weighted sum 



function we were looking at. 



□ 




■m 




3 = 1 
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where dxj =pk/qk - Pk+i/qk+i- 

Abel summation reduces this to the Birkhoff sum S m — J2T=i 9(^ a ) with g(x) = 
cot(7ra;) for which we have shown the limit to exist. 

The key is that for g(x) = cot(7rx) and derivatives, the Rolle error function is 
bounded. The monotonicity follows from the fact that /(a) + a 2 f is a csc 2 (x) 
Birkhoff sum, which is monotone. The limit of /„ is equal to the limit of / n (x) = 
Y/j=i x <ln\9(jln)dxj with 7„ = (p n+1 /q n+1 +p n /q n )/2. Let rj denote the Rolle 
points. Then 

n 

fn{x) =J2g(rj) dxj 
i=l 

Abel summation allows to replace this by 

fn(x) = S [xqn] — - -j S k ■ 

in q n k=1 

Because Sk < k, this converges. The first term is handled by the golden graph 
theorem. The second part is just a Cesaro average. 

The imaginary part of g is Hecke's Birkhoff sum g% (x) = 2ir{x} = 2tt(x~ [x] — 1/2). 
It leads to a limiting function satisfying f(ax) = a 2 f(x) because g" = 0. The Abel 
sum gives 

1 1 [X£n] 

fn(x) = [xq n ] 5- Y] 1 . 

In ql ^ 

This explains why we have exact selfsimilarity and smoothness of the limiting dif- 
ference. Also this was experimentally first noted in [11]. □ 

We could also modify the proof done for cot. Wc have again a function r(y) 
which is explicitly known. We write <j(y) — T(y) as a sum and show convergence. 
We can do the summation even so we have no definite sign because we know the t 
function and can use Abel summation to show that the limits exist. 
Why does secret;) = 1/ cos(7rx) give a different limiting golden graph than for cot? 
The reason is a different pole at x = 7r/2 and is discontinuous at 0. The trigono- 
metric identity just mentioned shows however that we can reduce it as a sum of 
two cot Birkhoff sums with different starting points. 



The function R(a) — J2ugz sec ( 7Tna ) satisfies R(l/a) — —aR(a) andi?(a+2) = 
R(a). It is a modular form of weight 1 and R — 9 2 , where 9(a) = 2^2 n ez e ™ n a 
V , y q"- 2 which is a classical 9 function. 1 To understand the curlicue problem, 
finite sums of the later would be important. 

As an application we are so able to deal with the Birkhoff sum studied in [TOl E] 
which was historically first studied in [B] who showed that 



1 " 

- V"log|csc(fea)| 



fc=i 



converges to log(2) for almost all a implying that Y^=i 1°S |2sin(fca)| = o(n). We 
have seen in [10] that this is 0(log(n) 2 ) for strongly Diophantine a and have seen 
that for Diophantine a of constant type, the Birkhoff sum is 0(log(n)). 
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Figure 9. The limiting Birkhoff sum function s'(x) 
linin-j.oo Se n with the Birkhoff sum 



f 



— y sec(nka) 



k=l 



has period 6. It is related to R — 9 2 with R(a) = 1+2 X^fcLi sec(a), 
where im(a) > 0. The attractor consists of 6 different limiting 
functions hin^^oo S6n+fc with k — 0, . . . , 5. We get the same at- 
tractor for cot(a; + n/2) = — tan(x) when starting from x = 0. 



The difference limit function we have seen in is the sum of a selfsimilar 
golden graph with different initial condition and a correction which comes from 
Abel summation and which is close to the graph of the esc Birkhoff limit of the s'. 
The limiting behavior seen in [TT] can be explained with the cot series. 

The following corollary is related to Lemmas 5 and 6 of Hardy-Littlewood's 1928 
paper [5], a paper which did not focus on the particular golden rotation numbers 
as we do. 



Corollary 3.4 (Cosecant Birkhoff sum) . For small \y\, the Birkhoff sums S\ xqn \{y) / = 

Y^k=i 9'{y/ln + ka) of g'(x) = — ir csc 2 (irx) converge to a monotone function s' y (x) 
along even and odd subsequences. For y = 0, the even and odd limits agree and the 
graph satisfies s'(ax) = as'(x). 



Proof. Differentiate the result for cot. 



□ 
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What class of functions g have the property that (l/q n ) J2k=i 9(ka) converges 
along some subsequence? We have seen that this is true for f{x) = c cot(irx) + h(x) 
where h is continuous. This includes virtually all functions which have a single 
pole at 0. A different type of pole appears with sec(7r:r) and also here, we have got 
limiting function, as seen in Figure (|9]). 

For some other Diophantine numbers like quadratic irrationals, the story is similar. 
Instead of a fixed point, we get a periodic attractor. In |13j . the expectation over 
a was taken, in which case one gets a limiting distribution and no subsequence has 
to be chosen. This suggests that for almost all a, the renormalized random walks 
converge to an attractor with a limiting distribution. For Liouville a, the situation 
is believed to be different and no limiting function s(x) can exist. Indeed we can 
construct a close to rational numbers so that the Birkhoff sum is unbounded. 
What about functions g with more than one pole? An example is g(x)+g(x+j3) with 
g(x) = cot(7ra). If j3 is rational, then a subsequence still converges and the graph is 
similar to a Birkhoff sum with one pole and quadratic irrational. For (3 rationally 
independent from 1 and a we do not know what happens. The situation is even 
more difficult then than for non Diophantine a, where we still have the continued 
fraction expansion. The challenge is to determine, how large the Birkhoff sum of 
g(na) + g(na + 0) can become. There is no reason why a Birkhoff limiting function 
should then exist. Indeed, experiments indicate that s m (x) explodes in general with 
a speed which depends on arithmetic properties of a and /3. 
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